UMN-TH-1328/95 
February 1995 

BIG BANG NUCLEOSYNTHESIS* 

KEITH A. OLIVE 

School of Physics and Astronomy, University of Minnesota 
Minneapolis, MN 55455, USA 
E-mail: 01ivc@mnhcp.hcp.umn.edu 



ABSTRACT 

The current status of big bang nucleosynthesis and its implications for physics 
beyond the standard model is reviewed. In particular, limits on the effective 
number of neutrino flavors and extra Z gauge boson masses are updated. 

1. Introduction 

The overall status of big bang nucleosynthesis is determined by the comparison 
of the rather slowly changing theoretical predictions of the light element abundances 
and the sometimes quickly changing observationally determined abundances. The 
observed elements, D, 3 He, 4 He, 7 Li, have abundances relative to hydrogen which 
span nearly nine orders of magnitude. By and large, these observations are consistent 
with the theoretical predictions and play a key role in determining the consistency of 
what we refer to as the standard big bang model and its extrapolation to time scales 
on the order of one second. Here, I will review the status of this consistency. I will 
begin by briefly outlining the key sequence of events in the early Universe which leads 
to the formation of the light elements. I will then discuss the current status of the 
observations in relation to theory of each of the light elements. Finally, I will discuss 
the current limits on physics beyond the standard model. 

2. A Brief Primer on the Theoretical Predictions 



Conditions for the synthesis of the light elements were attained in the early Uni- 
verse at temperatures T <, 1 MeV, corresponding to an age of about 1 second. At 
somewhat higher temperatures, weak interaction rates were in equilibrium, thus fix- 
ing the ratio of number densities of neutrons to protons. At T ^> 1 MeV, (n/p) ~ 1. 
As the temperature fell and approached the point where the weak interaction rates 
were no longer fast enough to maintain equilibrium, the neutron to proton ratio was 
given approximately by the Boltzmann factor, (n/p) ~ e _Am//T , where Am is the 
neutron-proton mass difference. The final abundance of 4 He is very sensitive to the 
(n/p) ratio. 
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The nucleosynthesis chain begins with the formation of deuterium. However, 
because the large number of photons relative to nucleons, r]~ l = n^jus ~ 10 10 , deu- 
terium production is delayed past the point where the temperature has fallen below 
the deuterium binding energy, Eb = 2.2 MeV. When the quantity 77 _1 exp(— Eb/T) ~ 
1, the nuclear chain begins at a temperature T ~ O.lMeV. 

The dominant product of big bang nucleosynthesis is 4 He resulting in an abun- 
dance of close to 25% by mass. In the standard model, the 4 He mass fraction depends 
primarily on the baryon to photon ratio, rj. The change due to the uncertainty in the 
neutron half-life is small (this effect is shown in Figure 1). When we go beyond the 
standard model, the 4 He abundance is very sensitive to changes in the expansion rate 
which can be related to the effective number of neutrino flavors as will be discussed 
below. Lesser amounts of the other light elements are produced: D and 3 He at the 
level of about 10~ 5 by number, and 7 Li at the level of 10 -10 by number. 



Fig. 1. The light element abundances from big bang nucleosynthesis. 

The resulting abundances of the light elements are shown in Figure 1 from the 
calculations in ref. 1. The curves for the 4 He mass fraction, Y, bracket the computed 
range based on the uncertainty of the neutron mean-life which has been taken as 
r n = 887 ± 2 s. The 4 He curves have been adjusted according to the corrections in 
ref. 2. Uncertainties in the produced 7 Li abundances have been adopted from the 
results in ref. 3. Uncertainties in D and 3 He production are negligible on the scale of 
this figure. The boxes correspond to the observed abundances and will be discussed 
below. It is clear that as the observational boxes line up on top of each other, there is 
an overall agreement between theory and observations in the range rjio = 10 10 r] = 2.8 
- 3.9. 

3. The Observations 

Because helium is produced in stars, it is very difficult to extract the primordial 
abundance from the observations. Ideally, one would want to look for primordial 
helium in regions where the stellar processing is minimal, i.e., in regions where the 
abundances of elements such as carbon, nitrogen and oxygen are very low. The 4 He 
abundance in very low metallicity regions is best determined from observations in 
extragalactic HII regions of Hell — > Hel recombination lines. There are extensive 
compilations of observed abundances of 4 He, N, and O, in many different galaxies Si. 
In Figure 2, the 4 He vs. O data is shown along with its associated linear fit i (details 
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Tabic 1. Linear Fits for Y vs. 0/H. 



Set 


# Regions 


r 


X 2 /dof 


Y P 


10 2 x slope 


1 P 


All 


49 


0.56 


0.78 


.234 ± 


.003 


1.14 ±0.24 


0.239 


1st cut 


41 


0.51 


0.61 


.232 ± 


.003 


1.38 ±0.36 


0.238 


-outliers 


34 


0.45 


0.70 


.232 ± 


.003 


1.39 ±0.38 


0.238 


2nd cut 


21 


0.41 


0.64 


.229 ± 


.005 


2.37 ± 1.13 


0.238 


-outliers 


19 


0.40 


0.70 


.229 ± 


.005 


2.42 ± 1.15 


0.238 


C 


22 


0.35 


0.71 


.232 ± 


.003 


1.58 ±0.54 


0.238 



of this fit are given in the second line of Table 1). 



Fig. 2. The observed abundances of 4 He vs. O/H in extragalactic HII regions along with a linear fit 
to the data. 

In Table 1, various fits to the data and subsets of the data are given. Details 
concerning the subsets of the data shown can be found in ref. 6. As one can see 
there is a considerable degree of stability in these fits, leading to a 2 a upper limit of 
0.238 for the primordial abundance of 4 He. There is in addition an overall systematic 
uncertainty of about 0.005 in Y p giving a range (2 a plus systematic) of 0.221 - 0.243 
for Y p and is shown in Fig. 1 as the large box bracketing the 4 He curves. 

It is more difficult to compare the primordial deuterium and 3 He abundances 
with the observations. Despite the fact that all observed deuterium is primordial, 
deuterium is destroyed in stars. A comparison between the predictions of the stan- 
dard model and observed solar and interstellar values of deuterium must be made 
in conjunction with models of galactic chemical evolution. The problem concerning 
3 He is even more difficult. Not only is primordial 3 He destroyed in stars but it is 
very likely that low mass stars are net producers of 3 He. Thus the comparison be- 
tween theory and observations is complicated not only by our lack of understanding 
regarding chemical evolution but also by the uncertainties of the production of 3 He 
in stars. 

It appears that D/H has decreased over the age of the galaxy. The (pre)-solar 
system abundance of deuterium is0 D/H rs (2.6 ± 1.0) x 10 -5 , corresponding to an 
age t ~ 9 Gyr, while the present (t ~ 14 Gyr) ISM abundance of D/H is! 1.65 x 10~ 5 . 
Thus, if 7710 is in the range 2.8 - 3.9, then the primordial abundance of D/H is 
between 4.5 - 8 xl0~ 5 , and it would appear that significant destruction of deuterium 
is necessary. Chemical evolution models which destroy D by as much as a factor of 5 
have been considered recently!. 
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Note that there is a reported detection of D in a high redshift, low metallicity 
quasar absorption system!^ with an abundance which may be the primordial one. 
This observation is shown in Fig. 1 by the small box on the D/H curve at a value 
of r/io m 1.5. As one can see the corresponding value of Y p (at the same value of 
77) is in excellent agreement with the data. 7 Li is also acceptable at this value as 
well. Due to the still some what preliminary status of this observation (in fact a 
recent reportS claims a much lower D abundance along a different line of sight) and 
the fact that it can also be interpreted as a H detection in which the absorber is 
displaced in velocity by 80 km s _1 with respect to the quasarS, it is premature to fix 
the primordial abundance with that value. A high value for the D abundance would 
require an even greater degree of D destruction over the age of the galaxy. 

There are however potential problems for 3 He. The lower limit on r\ was derivedll! 
by noting that although stars can destroy 3 He, even very massive stars, still preserve 
at least 25 % of the initial D + 3 He. (It is the sum of D and 3 He that is important as 
D is burned to 3 He in the premain-sequence phase of stars.) A value of r/10 lower than 
2.8 would yield (D+ 3 He)/H > 10~ 4 so that even if the maximal amount of 3 He is 
destroyed, it would still exceed the presolar valuei of (D+ 3 He)/H m 4.1 ± 1.0 x 10 -5 . 
But, in low mass stars, 3 He is produced rather than destroyed. A calculation!^ of the 
final 3 He abundance relative to the initial abundance of D + 3 He gives, 
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( 3 He/H) / = 1.8 x 10~ 4 + 0.7 [(D + 3 He)/H 

so that a 1 M star produces 2.7 times as much 3 He as the initial D + 3 He (for (D 
+ 3 He) initia i = 9 x 10~ 5 ). This would lead to an evolutionary behavioral of the type 
shown in Fig. 3. The chemical evolution model has been chosen so that D/H agrees 
with the data and assumes that t]\q = 3. The problem being emphasized concerns 
3 He and can be seen by comparing the solid curve with the filled diamonds. 



Fig. 3. The evolution of D/H (dashed curve), 3 He/H (solid curve) and (D + 3 He)/H (dotted curve) 
as a function of time. Also shown are the data at the solar epoch t « 9iLGyr and today for D/H 
(open squares), 3 He/H (filled diamonds) and (D + 3 He)/H (open circle) Ji3. 

A number of questions regarding the 3 He discrepancy can be raised. First, one can 
ask whether or not the formula for 3 He production is valid. There is indeed evidence 
that 3 He is produced in planetary nebulaeB where abundances are found as high as 
3 He/H ~ 10~ 3 . Secondly, one may ask how uniform are the 3 He measurements. As 
it turns out, they are in fact not very uniform^ and show variations as large as a 
factor of 5 between different HII regions. Indeed, there may even be a correlation 
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between the size (mass) of the region and the amount of 3 He observedEl. Thus it 
may be possible that the solar values are depleted in 3 He. It is hoped that future 3 He 
observations will help to resolve this puzzle. 

Finally I turn to 7 Li. Over the last several years, there has been a considerable 
increase in the number of 7 Li observations^. 7 Li in old, hot, population II stars, is 
found to have a very nearly uniform abundance. For stars with a surface temperature 
T > 5500 K and a metallicity than about l/20th solar, the 7 Li abundance shows little 
or no dispersion beyond what is consistent with errors of individual measurements. 
The corresponding mean 7 Li abundance is Li/H = 1.2 ± .1 x 10~ 10 . Systematic errors, 
however, dominate the uncertainty in the Li abundance, allowing a value Li/H <,2x 
10~ 10 . These values (and their uniformity) should be compared with observations of Li 
in younger stars where the abundance can be much larger (by an order of magnitude) 
and shows considerable dispersion. Two key questions remain, however: how much 
of the observed Li is primordial (since Li is known to be produced); and how much 
of the primordial Li remains in the stars where Li is observed? 

Aside from the big bang, Li is produced together with Be and B in cosmic ray 
spallation of C,N,0 by protons and a-particles. Li is also produced by a — a fusion. 
Be and B have recently been observed in these same pop II stars and in particular 
there are a dozen or so, stars in which both Be and 7 Li have observed. Thus Be 
(and B though there is still a paucity of data) can be used as a consistency check on 
primordial Li EI. Based on the Be abundance found in these stars, one can conclude 
that no more than 10-20 % of the 7 Li is due to cosmic ray nucleosynthesis leaving 
the remainder (an abundance near 10~ 10 ) as primordial. It is also possible however, 
that some of the initial Li in these stars has been depleted. Standard stellar models^ 
predict that any depletion of 7 Li would be accompanied by a very severe depletion of 
6 Li. Until recently, 6 Li had never been observed in hot pop II stars. The observation^ 
of 6 Li which turns out to be consistent with its origin in cosmic-ray nucleosynthesis 
and with a small amount of depletion as expected from standard stellar models is a 
good indication that 7 Li has not been destroyed in these starsil. 

Consistency of the standard model of big bang nucleosynthesis relies on the con- 
cordance between theory and observation of the light element abundances for a single 
value of rj. I now summarize the constraints on 77 from each of the light elements. 
From the 4 He mass fraction, Y < 0.238(0.243), we have that 7710 < 2.5(3.9) as a 2cr 
upper limit (the higher value takes into account possible systematic errors). Because 
of the sensitivity to the assumed upper limit on Y p , the upper limit on 77 from D/H, 
though weaker is still of value. From D/H > 1.5 x 10~ 5 we have r? 10 < 7. The lower 
limit on rj, comes from the upper limit on D + 3 He and is 7710 > 2.8 if one ignores 
3 He production. Finally, 7 Li allows a broad range for 77 consistent with other light 
elements. When both uncertainties in the reaction rates and systematic uncertainties 
in the observed abundances are taken into account, 7 Li allows values of 7710 between 
1.3 - 4.5. Taken all together, these bound on rj constrain the fraction of critical den- 
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sity in baryons to be between 0.01 - .1, for a hubble parameter, h , between 0.4 - 1.0 
(the corresponding range for VL B h 2 Q is 0.010 - 0.016). 

4. Constraints on Physics Beyond the Standard Model 

Limits on particle physics beyond the standard model are mostly sensitive to 
the bounds imposed on the 4 He abundance. As is well known, the 4 He abundance is 
predominantly determined by the neutron-to-proton ratio just prior to nucleosynthesis 
and is easily estimated assuming that all neutrons are incorporated into 4 He, 

2(n/p) 

P 1 + (n/p) {l) 

As discussed earlier, the neutron-to-proton ratio is fixed by its equilibrium value at 
the freeze-out of the weak interaction rates at a temperature Tf ~ 1 MeV modulo 
the occasional free neutron decay. Furthermore, freeze-out is determined by the 
competition between the weak interaction rates and the expansion rate of the Universe 



G/T f b ~ r wk (7» = H(T f ) ~ y GnNTj (3) 

where N counts the total (equivalent) number of relativistic particle species. The 
presence of additional neutrino flavors (or any other relativistic species) at the time 
of nucleosynthesis increases the overall energy density of the Universe and hence the 
expansion rate leading to a larger value of Tf, (n/p), and ultimately Y p . Because of 
the form of eq. (0) it is clear that just as one can place limit JhI on N, any changes 
in the weak or gravitational coupling constants can be similarly constrained (for a 
recent discussion see ref. 23. 

In the standard model, the number of particle species entering into eq. (|3]) can 
be written as N = 5.5 + \N U (5.5 accounts for photons and e ± ). The observationally 
derived primordial 4 He abundancJ of Y p = 0.232 ± 0.003 ± .005 translates into a best 
value for N u = 2.2±0.27±.42 which implies a 2a upper limit of 2.74 which is extended 
to N v < 3.16 when systematics are included. At face value, such a limit would 
exclude even a single additional scalar degree of freedom (which counts as such as 
a majoron unless it decoupled early enoug hi so that its temperature, Tb at the time 
of nucleosynthesis was suppressed and (T B /T U ) 4 < |(.16) = .28. In models with right- 
handed interactions, and three right-handed neutrinos, the constraint is more severe. 
The right-handed states must have decoupled early enough to ensure (T UR /T UL ) i < 
(.16)/3 ~ .05. The temperature of a decoupled state is easily determined from entropy 
conservation, (T x /T„) = ((A3/A)/N(T d )) 1 ^ 3 . One additional scalar requires N(T d ) > 
28 or decoupling must have occurred above the QCD phase transition at a temperature 
T d > Tq CD ~ 200 MeV. Three right-handed neutrinos would require N(T d ) > 100, 
so that T d > Mw- If right-handed interactions are mediated by additional gauge 
interactions, the associated mass scale becomes Mz> > O(100) TeV! 



6 



The limits on N u , however, are sensitive to the upper limit on 4 He which is in turn 
sensitive to assumed systematic errors and to the lower bound on rj. In addition, the 
limits described above may be overly restrictive^. The best value for N v is 2.2 and 
may in fact be unphysical if v T is lighter than ~ 1 MeV, as is quite likely. In this 
case, the limits on N v must be accordingly renormalizedEZl. In Fig. 4, the effect of 
renormalizing the limit on N v is shown. 



Fig. 4. The 95 % CL upper limit on N v as a function of the systematic uncertainty in Y p . The 
solid (dashed) curves correspond to the condition that N v > 2 (3) . Each of these two cases is shown 
for three choices of 7710: 1.5, 2.8, and 4.0. The smaller valueSpOf rj yield weaker upper limits and 
systematic errors have been described by a top-hat distributionEZl. 

In summary, I have argued for the overall agreement between theory and observa- 
tions as they pertain to the light element abundances as well as concordance between 
big bang nucleosynthesis and galactic cosmic-ray nucleosynthesis. There are however, 
open issues: Are the quasar line-of-sight measurements^ of D/H real; Why isn't there 
more 3 He, particularly in the solar system; Can the statistical and systematic errors 
in 4 He measurements be reduced; Can the large systematic errors in the 7 Li abun- 
dance be reduced? Clearly new data will be necessary to resolve these problems. 
Nevertheless, in spite of these uncertainties, nucleosynthesis is still able to set strong 
constraints on physics beyond the standard model. 

This work was supported in part by DOE grant DE-FG02-94ER40823. 
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